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ABSTRACT 


This  report  presents  a description  and  interpretation 
of  an  unusual  type  of  isolated  atmospheric  gravity  wave 
observed  near  Tennant  Creek  in  central  Australia.  Comparison 
of  experimental  data  with  theory  shows  that  waves  of  this 
type  belong  to  a new  class  of  deep  fluid  internal  solitary 
v/ave  discovered  by  Benjamin  and  by  Davis  and  Acrivos.  An 
interesting  feature  of  these  solitary  waves  is  that  they  mark 
the  onset  of  significant  atmospheric  turbulence.  It  is  proposed 
that  the  following  source  mechanisms  play  an  important  role  in 
the  creation  of  solitary  atmospheric  waves  in  the  arid  interior 
of  Australia  :(i)  the  impulsive  interaction  of  an  intense 
thunderstorm  downdraft  or  downdraft  generated  gravity  current 
with  an  inversion;  (ii)  the  interaction  of  nocturnal  katabatic 
flow  with  an  existing  radiation  inversion;  and  (iii)  the 
interaction  of  stationary  gravity  lee  waves  or  standing  eddies 
with  a developing  nocturnal  inversion  following  a reversal  in 
the  circulation  over  an  orographic  feature. 
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1.  INTRODUCTION 


It  is  well  known  that  the  creation,  propagation  and 
dissipation  of  gravity  waves  play  an  important  role  in  the 
dynamics  of  the  earth's  atmosphere.  The  reasons  for  studying 
long-period  waves  of  this  type  are  therefore  two-fold:  in  the 

first  place,  an  understanding  of  the  initial  disturbance  which  gives 
rise  to  these  waves,  the  mechanisms  by  which  they  propagate  and 
their  dispersion  characteristics  provides  insight  into  the 
fundamental  nature  of  atmospheric  fluid  mechanics;  secondly,  it 
seems  to  be  equally  important  to  understand  the  influence  that 
the  passage  of  these  waves  has  on  the  development  of  local 
meterological  conditions  - for  example,  gravity  waves  may  often 
be  correlated  with  the  release  of  latent  atmospheric  instabilities 
which  in  turn  give  rise  to  such  diverse  phenomena  as  clear  air 
turbulence,  propagating  squall  lines  and,  possibly,  in  the 
extreme  case  the  development  of  tornado  activity. 

The  purpose  of  this  report  is  to  summarize  the  principal 
results  of  a preliminary  study  of  subsonic  atmospheric  waves  which 
were  observed  during  the  period  from  August  1975  to  June  1976 
using  a five-component  array  of  high  sensitivity  microbarometers 
installed  at  the  Warramunga  Seismic  Station  located  near  Tennant 
Creek  in  central  Australia. 

Of  foremost  interest  to  this  study  are  the  observations 
of  large  amplitude  isolated  waves  of  permanent  form  which  are 
thought  to  represent  a new  type  of  naturally  occurring  internal 
solitary  wave.  This  report  will  be  mainly  concerned  with  a 
description  and  interpretation  of  these  unusual  waves.  A variety 


of  large  amplitude  long  period  dispersed  gravity  w ave  trains  have 
also  been  observed  during  this  period.  Discussion  of  these  events 
will  be  deferred  until  a later  report. 

2.  DATA  RECORDING  AND  PROCESSING  TECHNIQUES 

The  Warramunga  Seismic  Station  is  situated  37  kilometers 
SSE  of  Tennant  Creek.  The  centered  quadrilateral  array  of 
microbarometers  is  located  on  slowly  undulating  semi-desert 
terrain  at  an  elevation  of  about  4l0  meters  MSL  with  relief  in 
the  area  encompassed  by  the  array  rising  to  a maximum  of  about 
6 meters.  The  nearest  significant  landscape  feature  is  a 
barren  50  meter  granite  ridge  which  runs  for  about  200  meters  in 
a north-south  direction  about  1 kilometer  east  of  site  1. 

Probably  the  most  important  topological  features,  insofar 
as  the  location  of  this  array  is  concerned,  are  the  600  meter 
Murchison  and  Davenport  Ranges  which  run  from  20  to  160  kilometers 
to  the  SSE  and  the  1600  meter  Macdonnell  Ranges,  450  kilometers 
to  the  south.  All  of  the  land  to  the  west  of  the  array  can  be 
described  as  semi-featureless  stoney  desert  and  an  extensive 
dry  steppe  area  known  as  the  Barkly  Tablelands  forms  the  north- 
east quadrant.  Further  afield,  the  major  orographical  features 
are  the  5000  meter  mountains  of  central  New  Guinea,  1800  kilometers 
to  the  NNE,  and  the  2400  meter  Australian  Alps  of  Victoria  and 
New  South  Wales,  2300  kilometers  to  the  SE. 

The  configuration  of  the  five-component  infrasonic  array 
along  with  the  location  of  a vault  containing  a vertical  long- 
period  seismometer  are  illustrated  in  Figure  1.  An  array  of 
Daniels  type  noise  reducing  space  filters  arranged  in  the  form 
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of  a cross  v/ith  the  inlet  port  at  the  center  has  been  installed 
at  each  site.  These  filters  provide  useful  suppression  of 
incoherent  wind  noise  in  the  period  range  below  20  seconds. 

A least-squares  fit  to  the  measured  amplitude  response  of  the 
microbarometers  is  shown  in  Figure  2. 

Recently,  a SIAP  2000  weather  station  which  provides  a 
continuous  record  of  wind  speed  and  direction,  temperature, 
relative  humidity  and  precipitation  was  installed  midway 
between  the  microbarometer  at  site  1 and  the  long-period 
seismometer  vault. 

The  five  infrasonic  channels  and  the  vertical  long-period 
seismometer  channel  are  sampled  at  a rate  of  2 samples  per 
second,  digitized  via  a 14-bit  ADC  and  recorded  in  IBM 
compatible  format  on  7-track,  556  bpi  tape.  At  the  Australian 
National  University  these  tapes  are  converted  to  9-track,  800 
bpi  in  a compact,  2 data  word  per  24-bit  computer  word,  format 
which  conserves  magnetic  tape  and  which  is  suitable  for  analysis 
on  a 48K  Harris  Datacraft  6024/4  computer. 

The  main  signal  processing  technique  used  on  the  array 
data  consists  of  a beam-forming  program  which  utilizes  the 
non-linear  N-root  method  devised  by  Muirhead  (1968).  The 
output  ei(t)  of  the  i-th  array  element  is  digitally  filtered 
in  the  bandpass  region  of  interest,  properly  phased  by  time 
shifting  to  correspond  to  the  passage  of  a plane  wavefront 
across  the  array,  and  then  reduced  to  the  N-th  root  with  sign 
preserved.  The  resulting  product  is  summed  over  the  n-element 
array  to  give 


R„(t) 


n 


£ 
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ei(t ) 


1 

N 


signum  ( ei (t ) ) 


1 
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(2.1) 
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w’nere  e^(t)  represents  the  filter  output.  This  quantity  is 
then  raised  to  the  N-th  power  with  sign  preserved  to  provide 
the  signal  statistic 

N 

SN(t)  = | Rn( t ) | . signum  (R^Ct))  (2.2) 

The  principal  advantage  of  this  type  of  automatic 
infrasonic  array  processing  arises  from  the  fact  that  this 
non-linear  beam-forming  technique  gives  added  weight  to  the 
presence  of  coherent  energy  in  the  spectrum.  This  has  proven 
to  be  of  value  in  the  treatment  of  data  from  Tennant  Creek 
due  to  the  prevalence  in  this  semi-desert  environment  of 
incoherent,  essentially  non-Gaussian  noise  such  as  that 
caused  by  very  small  whirl-winds  (dust  devils)  which  interact 
with  only  one  sensor.  It  can  be  easily  shown  that  large 

events  of  this  type  are  suppressed  in  the  output  of  the  N-th 

-N 

root  process  by  a factor  of  approximately  n 

All  infrasonic  data  is  processed  for  a value  of  N=2  in 
order  to  take  advantage  of  a hardware  square  root  SAU  feature 
of  the  Datacraft  computer.  This  output  is  supplemented  by  a 
parallel  determination  of  an  integrated  polarity  stack  over 
the  phased  array  which  corresponds  to  the  N-th  root  process 
in  the  limit  N-*».  Further  details  on  the  N-th  root  multi- 
channel filter  may  be  found  in  Kanasewich  et  al.  (1973)  and 
Muirhead  and  Ram  Datt  (1976). 

Power  spectral  estimates  of  the  data  are  computed  directly 
using  the  Fast  Fourier  Transform  (Cooley  and  Tukey , 1965). 

In  order  to  reduce  the  variance  of  the  estimate  the  calculations 
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are  carrier!  out  using  ths  method  of  time  averaging  over  modified 
periodograms  described  by  Welch  (1967).  The  digital  time 
series  t is  divided  into  K portions  each  of  length  M=210 
points;  the  elements  in  each  segment  i are  then  weighted 
according  to  a cosine  taper  data  window  d^  effective  over 
the  10*  limits  of  the  segment,  as  suggested  by  Bingham  e_t  al . 
(1967),  and  transformed  to  the  frequency  domain  to  give 


xl(r)  - E 

— n * n j 


(2.3) 


The  power  spectral  density  is  then  determined  by  averaging 
over  the  K segments. 


TU  I Xi^f^  I » 


(2.4) 


where  At  is  the  sampling  interval  and  the  factor 


1 11  2 

0 M ET  ' °-875 


(2.5) 


preserves  the  invariance  of  the  area  under  the  spectrum  to 
the  influence  of  the  data  taper  window.  The  power  spectral 
estimates  obtained  in  this  way  are  converted  to  true  spectral 
estimates  by  dividing  by  the  square  of  the  microbarograph 
amplitude  response  shown  in  Figure  2. 
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3-  EXPERIMENTAL  RESULTS  AND  INTERPRETATION 
3- I Internal  Solitary  Atmospheric  Waves 

The  phenomenon  of  the  solitary  wave,  by  its  essentially 
nonlinear  nature,  occupies  a unique  place  in  the  development 
of  the  theory  of  wave  propagation  in  fluids.  These  waves  are 
usually  defined  (Lamb,  1932,  §252)  as  waves  of  single  elevation 
which  propagate  at  uniform  velocity  without  change  of  form. 

They  arise  from  a balance  between  the  tendency  of  waves  to 
steepen  ahead  of  their  crests  (amplitude  dispersion)  and  the 
tendency  of  the  Longer  period  wave  components  to  propagate  at 
higher  velocities  (frequency  dispersion). 

The  classical  solitary  wave  of  elevation  was  first  observed 
by  Russell  (1844)  on  the  free  surface  of  shallow  water  of  uniform 
depth.  Subsequently,  Boussinesq  (1871)  and  Rayleigh  (1876) 
independently  derived  approximations  for  the  speed  of  propagation, 
c,  and  the  form,  n(n),  of  the  wave  profile  at  the  free  surface. 
These  authors  found,  to  first  order  in  a = a/h. 


c2  = gh( 1+a) 


(3  .la) 


and 


(3-lb) 


where  h is  the  fluid  depth,  a is  the  maximum  amplitude 
of  the  v/ave  and  g is  the  acceleration  due  to  gravity.  These 
results  reveal  an  important  property  of  solitary  waves;  that 
is,  they  are  always  supercritical  - the  speed  of  propagation 


walways  exceeds  the  speed 


- 1 


of  infinitesimal  long  waves.  This  property  may  be  expressed 
to  first  order  in  the  dimensionless  amplitude  as 

P2  = 1 + a > 1 (3.3) 

where  F is  the  Froude  number.  It  should  be  noted  that  the 
classical  solitary  wave  is  derived  on  the  simplifying  assumption 
that  the  horizontal  length  scale  of  the  notion  is  long  compared 
to  the  fluid  depth. 

Korteweg  and  DeVries  (1895)  were  able  to  demonstrate,  in 
the  first  approximation,  the  existence  of  a class  of  periodic 
long  waves  of  finite  amplitude  and  permanent  form  which  are 
described  by  the  Jacobi  elliptic  function  aCn2(x/X,k)  of 
modulus  k.  In  particular,  they  showed  that  these  waves,  which 
they  called  onoidal  waves,  reduce  to  the  classical  solitary  wave  of 
finite  extent  in  the  limit  as  the  wavelength,  X , goes  to 
infinity.  Despite  this  progress,  it  was  not  until  relatively 
recently  that  Lavrentiev  (1946)  and  Friedrichs  and  Hyers 
(1954)  succeeded  in  rigorously  proving  the  mathematical 
existence  of  a solitary  wave  solution  of  the  full  non-linear 
equations . 

There  have  been  a number  of  attempts  to  improve  on  the 
first  order  solution  for  the  classical  solidary  wave  (e.g.. 

Long  (1956b) ;Laitone  (i960);  Byatt-Smith  (1970);  Grimshaw 
(1971);  and  Fenton  (1972)).  The  most  accurate  results  are  those 
of  Fenton,  who  obtained  a ninth-order  solution,  and  Byatt-Smith 
who  numerically  solved  for  the  case  of  the  solitary  wave  an 


exact  integro-differential  equation  for  surface  waves  of 
permanent  form.  The  principal  result  of  these  investigations 
is  that  the  solitary  wave  of  maximum  amplitude  corresponds 
to 


a = 0. 85  , F 
max  J 3 max 


= 1.31  - 


(3-4) 


It  is  worth  noting  that  the  first  approximation  to  the  shape 
of  the  classical  solitary  wave  (3. lb)  compares  favourably 
with  the  exact  profile  computed  by  Byatt-Smith  provided  the 
non-dimensional  amplitude  a is  less  than  about  0.7. 

The  theoretical  study  of  internal  solitary  waves,  that 


is,  waves  which  exist  as  a consequence  of  internal  density 
stratification,  was  initiated  by  Keulegan  (1953)  who  considered 
a system  of  two  superimposed  liquids  of  different  constant 
densities,  bounded  above  and  below  by  rigid  surfaces.  Long 
(1956b)  has  also  investigated  this  particular  model.  Abdullah 
(1956)  obtained,  using  a systematic  perturbation  procedure 
developed  by  Friedrichs  (1948)  and  Keller  (1948),  a first 
order  solution  for  the  classical  internal  solitary  wave  at 
the  interface  of  a two-layer  atmosphere  subject  to  the  condition 
that  the  hydrostatic  lav;  holds  for  the  upper  layer,  that  is, 
that  motion  of  the  free  surface  can  be  neglected.  This  solution 
is  a supercritical  wave  of  elevation  described  by 


n(x) 


a sech2 


(3-5a) 


and 

with 


c = cq ' ( l+^a ) 

4 log  (l+v/2) 
v;a,  = h 

3 n. 


(3.5b) 

(3.5c) 
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w’nere  the  infinitesimal  internal  wave  critical  speed  is  given 
by 


with 


c0'  = Vi7^ 

g(Pl~P2 ) 

g'  = 

Pi 


(3.6) 


(3.7) 


Here,  p2  and  pj  are  the  densities  of  the  upper  and  lower 
fluids  respectively,  h is  the  depth  of  the  undisturbed  lower 
layer  and  wa  is  the  full  width  of  the  wave  profile  at  half 
maximum. 

A more  general  treatment  of  internal  solitary  waves 
has  been  given  by  Peters  and  Stoker  (I960).  They  examined  the 
full  problem  of  a two-fluid  system  with  a free  upper  boundary 
as  well  as  the  more  difficult  problem  of  the  existence  of 
internal  solitary  waves  in  a fluid  whose  density  decreases 
exponentially  with  height.  In  the  case  of  the  two  fluid 
problem  they  found  two  solitary-wave  solutions;  one  type, 
with  maximum  amplitude  at  the  free  surface,  corresponds  to 
the  ordinary  classical  solitary  wave  of  elevation;  the  other 
type  describes  an  internal  solitary  wave  with  a maximum  amplitude 
at  the  interface  much  larger  than  the  amplitude  at  the  surface. 
The  internal  solitary  wave  may  be  either  a wave  of  depression 
or  a wave  of  elevation.  If  the  depth  of  the  undisturbed  lower 
fluid,  h,  is  less  than  half  the  depth  of  the  upper  fluid, 
the  streamlines  of  the  disturbance  are  lines  of  elevation  at 
the  interface  and  depression  at  the  free  surface;  if,  in  addition 
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the  difference  in  the  two  fluid  densities  is  small,  the 
solution  then  corresponds  to  the  approximate  form  given  by 
(3-5)  above.  This  approximation  will  be  assumed  to  provide 
a sufficiently  accurate  description  of  classical  internal  solitary 
waves  in  the  atmosphere  for  the  present  purposes. 

In  the  case  of  a fluid  of  finite  depth  with  a density 
distribution  which  decreases  exponentially  upwards , Peters  and 
Stoker  found  an  infinite  number  of  possible  internal  solitary 
wave  modes  corresponding  to  an  infinite  spectrum  of  internal 
critical  speeds.  This  problem  has  also  been  Investigated  by 
Long  (1965)  and,  in  an  elegant  and  very  general  treatment  of 
the  subject  of  internal  solitary  waves  in  shallow  fluids,  by 
Benjamin  (1966).  The  general  solitary  wave  solutions  in  this 
case  are  complicated  and  will  not  be  given  here;  it  is  noted, 
however,  that  according  to  Benjamin  the  solitary  wave  mode 
corresponding  to  the  highest  critical  speed  - predominantly  a 
wave  of  depression  - is  probably  the  most  significant  mode  in 
fluids  of  this  type. 

Up  to  this  point  this  discussion  has  been  limited  to  a 
consideration  of  internal  solitary  waves  which  can  exist  in 
fluids  of  finite  depth.  Benjamin  (1967)  and,  independently, 

Davis  and  Acrivos  (1967)  have  presented  the  results  of  a 
theoretical  and  experimental  investigation  of  an  entirely  new 
class  of  internal  solitary  wave  phenomenon  which  can  exist  in 
fluids  of  great  depth.  These  new  types  of  internal  waves  are 
shown  to  exist  in  regions  where  the  fluid  density  varies  only 
within  a layer  of  thickness  h which  is  much  smaller  than  the 
total  fluid  depth  and  smaller  than  the  effective  horizontal 
length  scale,  X , characteristic  of  the  solitary  wave. 
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For  a tv/o  fluid  system  in  which  the  upper  fluid 
extends  to  infinity  over  a lower  fluid  of  depth  h resting 
an  a horizontal  rigid  surface,  Benjamin  finds  a solitary  wave 
solution  of  the  form 

aA  2 

n(x)  = (3-8a) 

x2  + X2 


with 


where 


8 ^ h 
3 P2  a 


(3.8b) 

(3.8c) 


and  p2  and  p2  are  the  densities  of  the  upper  and  lower  fluids 
respectively. 

Benjamin  has  also  considered  the  problem  of  a shallow 
fluid  with  exponential  density  gradient  in  contact  with  a 
deep  fluid  of  constant  density,  and,  further,  the  more  general 
problem  of  solitary  waves  associated  with  a thin  transition 
region  contained  between  two  homogeneous  fluids  of  substantial 
depth.  In  direct  analogy  with  the  shallow  fluid  theory,  the 
continuous  variation  in  density  is  shown  to  give  rise  to  an 
infinite  set  of  solitary-wave  modes.  For  the  case  of  a thin 
region  of  sudden  vertical  density  variation  in  a stably 
stratified  system,  the  lowest-order  mode  solitary  wave  solution 
takes  the  form  of  a bulge  which  propagates  along  the  interface. 

This  interesting  result  has  been  confirmed  in  the  laboratory 
experiments  of  Davis  and  Acrivos  (1967)  and  Hurdis  and  Pao  (1975). 
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Consider  now  the  existence  of  solitary-wave  motions 

in  the  atmosphere.  It  would  seem  that  the  deep  fluid  solitary 
wave  solution  (3-8)  found  by  Benjamin  (1987)  is  best  suited 
to  a simple  description  of  solitary  waves  in  the  planetary 
boundary  layer  'while  classical  solitary  wave  theory  (3-5) 
should  provide  a reasonable  description  of  higher  altitude 
waves.  The  possible  existence  of  atmospheric  internal 
solitary  waves  was  suggested  by  Abdullah  (19^9). 

In  order  to  clarify  the  following  discussion  consider 
the  properties  of  internal  solitary  waves  associated  with  a 
density  discontinuity  at  the  950  mbar  (580m)  level  corrresponding 
to  a temperature  inversion  of  5°K,  as  predicted  by  both  the 
shallow  and  deep  fluid  theories.  The  air  beneath  the  inversion 
is  taken  to  have  a mean  temperature  of  270°K  and  the  dimensionless 
amplitude  of  the  wave,  a , is  chosen  to  have  a value  of  0.5, 
corresponding  to  a wave  with  amplitude  less  than  maximum.  With 
these  conditions,  a classical  internal  solitary  wave  as 
described  by  (3>5)  would  propagate  along  the  inversion 
at  a speed  of  12.6  m/sec  with  a full  width  at  half  maximum 
of  1.61  km  and  would  produce  a maximum  pressure  perturbation, 

AP,  at  ground  level  of  0.615  mbar;  alternatively,  for  the  same 
amplitude,  Benjamin's  theory  (3.8)  describes  a solitary  wave 
with  Wj  = 3-03  km  orooagating  at  a speed  of  11.8  m/sec.  At  the 
980  mbar  level  (280  m) , for  an  inversion  of  the  same  intensity, 
a similar  solitary  wave  would  produce  a ground  level  perturbation 
of  O.317  mbar  and  would  propagate  at  a velocity  of  8. 91  m/sec 

with  w,  = 0.806  km  according  to  the  classical  theory,  or  at  a 

'2. 

velocity  of  8.36  m/sec  with  w ^ = 1.52  km  as  given  by  the  deep 
fluid  theory. 
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A comparison  of  the  shallow  and  deep  fluid  solitary 
wave  profiles  for  an  Inversion  at  the  950  mbar  level  is  shown 
in  Figure  3-  Note  that,  as  a consequence  of  the  short  range 
exponential  rise  and  decay,  the  influence  of  the  classical  wave 
is  much  more  localized  than  that  of  the  deep  fluid  wave  of 
elevation.  The  shape  of  both  types  of  solitary  wave  depends 
on  the  value  of  ohe  dimensionless  amplitude  a ; as  a increases 
the  wave  profile  narrows  and  the  curvature  at  the  crest 
increases  until  at  maximum  amplitude  the  wave  is  reduced  to  a 
sharp  peak  enclosing  an  angle  of  120°  (Stokes,  1880;  Byatt- 
Smith,  1970).  This  variation  in  profile  is  illustrated  for  values 
of  a up  to  0.7  by  the  two  families  of  curves  shown  in  Figure  4 
with  the  depth  of  the  undisturbed  fluid  taken  as  the  unit  of 
length . 

3 . 2 Observations 

A large  number  of  isolated  waves  have  been  observed  at 
the  infrasonic  array  near  Tennant  Creek.  These  unique  waves, 
which  only  occur  at  night,  are  reasonably  well  described  by 
the  solitary-wave  theory  outlined  above.  Consequently,  as 
mentioned  in  the  introduction,  these  waves  are  interpreted 
as  being  examples  of  atmospheric  internal  solitary  waves 
which  propagate  along  the  noctural  inversion. 

Gossard  and  Hunk  (1954),  in  a study  of  atmospheric 
pressure  oscillations,  noted  several  examples  of  a gravity 
wave  train  preceeded  by  a positive  pressure  pulse.  These  early 
measurements  suffer  from  a lack  of  resolution  and  as  such  are 
not  directly  compared]  o with  the  observations  made  at  Tennant 
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Creek.  Nonetheless,  on  the  basis  of  available  evidence,  it 
seems  reasonable  to  assume  that  these  pressure  pulses  are  of 
the  same  type  as  the  isolated  waves  described  here. 

McAllister  et  al.  (1969)  have  reported,  in  a pioneering 
atmospheric  acoustic  sounding  study,  observations  of  "weak, 
front-like"  nocturnal  disturbances  characterized  by  a sharp 
spike  in  the  sounder  record  followed  immediately  by  a rapid 
buildup  of  turbulence  and,  after  a period  of  about  10  minutes, 
by  well-developed  sinusoidal  motion  in  slowly  separating 
strata.  These  observations  were  carried  out  during  the  month 
of  June  at  Ivy  Tanks  on  the  edge  of  the  arid  Nullarbor  Plain  in 
South  Australia.  These  authors  interpreted  these  spikes  as 
reflections  from  the  turbulent  core  of  a somewhat  broader 
disturbance . 

A very  interesting  series  of  acoustic  sounding  experiments 
carried  out  again  during  the  month  of  June  at  Julia  Creek, 
Queensland  have  been  described  by  Reynolds  and  Gething  (1970). 

On  several  occasions  they  also  observed  transient  sharp  spikes 
on  the  interface  of  the  nocturnal  inversion.  Some,  but  not  all 
of  these  spikes  were  associated  with  a jump  in  the  height  of 
the  inversion.  These  authors  noted  that  waves  associated  with 
these  sharp  spikes  on  the  inversion  level  might  represent  a 
form  of  non-linear,  large  amplitude  solitary  wave. 

There  is  little  doubt  that  the  solitary  waves  observed 
at  the  Tennant  Creek  infrasonic  array  are  of  the  same  type  as 
those  observed  800  kilometers  to  the  east  at  Julia  Creek  and 
1300  kilometers  to  the  south  at  Ivy  Tanks.  It  does  not  appear 
likely  at  this. point  that  these  waves  preserve  their  identity 
over  these  distances;  however,  this  possibility  should  be 


borne  in  mind. 
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Smart  (1966)  and  Jordan  (1972)  have  reported  observations 
of  "exponential  pressure  pulses"  near  Denver,  Colorado,  which 
are  similar  in  many  respects  to  some  of  the  solitary  waves 
described  in  this  report. 

The  only  other  evidence  for  atmospheric  solitary  waves 
appears  to  be  the  description  by  Abdullah  (1955)  of  a large 
amplitude  propagating  disturbance  which  appeared  over  Kansas 
during  the  early  daylight  hours  of  June  29,  1951.  This 
disturbance,  which  produced  a ground  level  pressure  perturbation 
of  3-4  mbar,  took  the  form  of  an  elevated  mass  of  cold  air 
propagating  on  an  inversion  at  a height  of  about  2 kilometers. 
The  elevated  disturbance  appeared  to  extend  over  about  150 
kilometers  and  was  observed  to  travel  with  approximately 
constant  form  at  speeds  between  18  m/sec  and  2|nr,/sec  over  a 
distance  of  about  300  kilometers.  Abdullah  concluded  that  this 
disturbance  represented  an  internal  solitary  wave  and 
attributed  its  formation  to  the  impulsive  movement  of  a 
quasistationary  cold  front  into  the  layer  of  inversion. 

As  can  be  seen  from  the  following  description,  the 
atmospheric  solitary  waves  observed  at  Tennant  Creek  appear  to 
be  entirely  different  from  the  single  observation  described  by 
Abdullah.  These  subsonic  waves  are  characterized  by  an  initial 
slow  increase  in  pressure  over  a period  of  several  minutes  which 
gradually  develops  into  an  exponential  increase  leading  to  a 
fairly  sharp  crest  at  which  point  the  pressure  falls  rapidly  to 
the  ambient  level.  A wide  variety  of  these  commonly  occurring 
gravity  waves  are  illustrated  in  Figures  5 to  10. 

In  most  cases,  these  waves  occur  as  isolated  waves  of 
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elevation  as  shown  in  Figures  5 and  6.  They  are  usually 
asymmetrical  but  a few  examples  have  been  observed  which  are 
almost  symmetrical  in  form.  Occasionally,  the  initial  positive 
pressure  pulse  Is  followed  by  a wave  of  depression  as  shown  by 
the  example  in  Figure  7(a).  A number  of  events  have  been 

observed  (see  Figures  7 to  9 ) in  which  the  initial  positive 
pressure  pulse  is  followed  by  several  long  period  large  amplitude 
components  which  may  be  either  coherent  or  incoherent  over  the 
array  - these  waveforms  are  similar  to  those  observed  by 
Smart  (1966)  and  Jordan  (1972).  Since  all  of  these  atmospheric 
pressure  oscillations  are  initiated  by  a positive  pressure  pulse 
of  the  same  form,  these  events  are  all  believed  to  be  manifest- 
ations of  the  same  phenomena. 

The  salient  features  of  these  atmospheric  solitary  waves 
may  be  seen  from  an  examination  of  the  summary  of  characteristics 
presented  in  Table  I. 

Waves  of  this  type  have  been  observed  with  amplitude  as 
high  as  700  dynes/cm2.  They  propagate  at  speeds  between  4 
and  lo  m/sec,  with  an  effective  wavelength,  as  measured  by  the 
full  width  at  half  maximum,  between  0.11*  kilometers  and  2.68 
kilometers  and  tend  to  occur  on  successive  nights  in  groups  of 
up  to  6 days  duration  as  shown  in  Figure  11.  Even  though  the 
available  data  spans  a period  of  only  10  months,  it  seems  clear 
that  the  frequency  of  occurrence  of  these  waves  is  seasonal. 

In  contrast  to  the  observations  reported  by  Smart  and  by  Jordan, 
they  do  not  occur  during  the  daylight  hours  between  08.30  and 
17.00  C.S.T.  (see  Figure  12).  A plot  of  the  frequency  of 
of  observation  as  a function  of  azimuth  (Figure  13)  reveals  the 
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fact  that  those  solitary  waves  originate  in  the  directions  of 
2C°  and  140  as  measured  from  true  north. 

An  interesting  feature  of  these  solitary  waves  is  their 
apparent  influence  on  ambient  atmospheric  conditions.  This 
feature  is  illustrated  in  the  diurnal  variations  of  the  micro- 
pressure spectrum  shown  in  Figure  14.  The  last  diagram  in 
this  series  (c)  has  been  chosen  to  illustrate  a typical 
40  hour  period  during  which  gravity  wave  activity  was  absent. 

As  can  be  seen,  the  late  evening  and  early  morning  hours  are 
characterized  by  a relatively  low-level  incoherent  long-period 
noise  spectrum  indicative  of  a quiet  stably  stratified 
noctural  atmosphere;  in  contrast,  the  onset  of  daybreak  marks  a 
transition  to  a state  of  turbulence  as  the  nocturnal  inversion 
is  dissipated  by  the  growth  of  a convectively  unstable  penetrating 
layer  at  the  surface.  This  transition  is  also  accompanied  by 
the  onsec  of  a light  surface  wind  of  about  4 m/sec. 

Consider  now  the  influence  of  the  presence  of  solitary 
waves  propagating  on  the  nocturnal  inversion  as  is  shown  in  the 
first  two  diagrams  in  Figure  14.  In  each  case,  the  passage  of 
the  isolated  wave  marks  the  onset  of  turbulence  in  the  atmosphere. 
This  remarkable  transition  is  particularly  evident  following  the 
passage  of  the  second  event  in  Figure  14(b)  and  may  be  seen  in 
greater  detail  in  Figure  6(a). 

A change  in  the  atmospheric  spectral  character  following 

/ 

the  passage  of  a solitary  wave  is  a feature  which  is  common  to 
all  events  of  this  type  which  have  been  observed  to  date  in 
Tennant  Creek;  however,  the  degree  to  which  this  phenomenon 
occurs  varies  over  wide  limits  as  can  be  seen  from  the  comparisons 
shown  in  Figure  15  of  power  spectral  densities  evaluated  before 
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and  after  the  passage  of  selected  events.  The  spectra 
computed  from  data  recorded  prior  to  the  passage  of  the  solitary 
wave  are  all  quite  similar  and  are  characteristic  of  low  noise 
stable  conditions  as  evidenced  by  the  presence  of  a ^ to  8 second 
period  peak  in  the  spectrum  due  to  nicrobarom  activity.  In 
examples  (a)  and  (b),  the  solitary  wave  marks  the  introduction 
of  spectral  components  over  the  complete  period  range  - it  is 
noteworthy  that  the  power  from  10  to  100  seconds  in  example  (a) 
is  increased  by  nearly  4 orders  of  magnitude.  In  examples  (c) 
and  (d)  additional  high  frequency  components  are  almost  completely 
absent  after  the  passage  of  the  wave;  example  (d)  is  interesting 
in  that  the  spectral  characteristics  before  and  after  the 
event  are  almost  identical  except  for  a small  contribution  in 
the  10  to  100  second  band. 

The  problem  posed  by  the  asymmetrical  form  of  these  isolated 
pressure  pulses  will  now  be  considered.  At  first  glance,  it 
would  seem,  on  the  basis  of  the  theory  reviewed  above,  that  if 
these  pressure  distributions  represent  an  internal  solitary  wave 
of  elevation,  then  they  should  be  symmetrical.  However,  all 
of  the  theoretical  treatments  are  limited  to  simple  static 
models.  Perhaps  the  most  serious  criticism  of  these  models 
when  applied  to  the  earth's  atmosphere  is  that  the  influence  of 
wind  shear  across  the  interface  (within  the  limits  of  dynamic 
stability)  is  neglected;  this  might  be  the  reason  for  the  observed 
asymmetrical  shape.  It  is  also  possible  that  the  turbulent 
wake  and  form  of  these  pressure  pulses  result  from  isolated  waves 
breaking  in  the  backwards  direction  - this  process  is  enhanced  by 
the  presence  of  a suitable  shearing  layer  (Long,  1956a).  This 
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explanation  does  not  appear  to  be  likely  but  it  should  be  borne 
in  mind. 

Another  explanation  of  these  wave  forms  is  indicated  by 
the  experimental  observations  and  numerical  calculations  of  Davis 
and  Acrivos  (1967)  on  solitary  v/aves  propagating  in  a thin 
layer  of  fluid  contained  between  two  deep  fluids  comprising  a 
stably  stratified  system.  They  found  that  waves  of  large  amplitude 
develop  closed  streamlines  characteristic  of  the  circulation  in  a 
vortex  pair  and  further,  that  these  waves  shed  semiperiodic 
v/aves  behind  the  main  disturbance.  It  follows  directly  from 
an  argument  based  on  dynamic  symmetry  (Benjamin,  1967)  that  a 
single  vortex  occurs  within  large  amplitude  solitary  waves  associated 
with  an  inhomogeneous  shallow  layer  of  fluid  lying  beneath  a 
much  deeper  homogeneous  fluid  and  that  this  internal  circulation 
could  lead  to  the  development  of  a turbulent  wake  as  is  observed. 

A final  possibility  which  could  account  for  the  observed 
wave  shape  is  that  these  large  amplitude  waves  significantly 
displace  the  fluid  above  the  inversion  level  and  thus  they  may 
effect  the  release  of  any  latent  atmospheric  instability.  The 
resulting  convective  mixing  would  last  for  a considerable  period 
of  time,  as  observed,  and  could  be  expected  to  result  in  a 
rapid  drop  in  pressure  to  the  ambient  level.  Some  support  for 
this  idea  is  provided  by  a comparison  (Figure  1 6)  of  the  pressure 
distribution  associated  with  a wave  which  initiated  only  slight 
changes  in  the  atmospheric  spectral  character  with  the  more 
typical  distribution  associated  with  a wave  with  a significantly 
turbulent  wake.  The  profile  of  the  wave  with  the  quiet  wake  is 
clearly  more  symmetrical.  Hov/ever,  this  interpretation  is 


-20- 


opan  to  criticism  since  some  of  the  observed  isolated  waves 
with  quiet  wakes  are  strongly  asymmetric. 

It  is  worth  noting  at  this  point  that  the  form  of  the 
leading  edge  of  these  solitary  waves  indicates,  as  could  be 
expected  from  the  atmospheric  scale  involved,  that  these  waves 
are  best  described  by  the  deep  fluid  internal  solitary  wave 
theory.  It  seems  that  a final  interpretation  of  these  unique 
waves  which  includes  an  explanation  of  the  asymmetric  profile 
will  involve  more  than  one  of  the  effects  discussed  above. 

It  must  be  emphasized  at  this  point  that  any  discussion 
of  the  source  mechanisms  which  generate  these  solitary  waves 
is  largely  conjectural  since  the  range  and  lateral  extent  of 
these  waves  are  almost  totally  unknown.  Existing  evidence,  such 
as  the  fact  that  these  waves  traverse  the  4 kilometer  infrasonic 
array  without  noticeable  attentuation  as  highly  coherent  planar 
wavefronts,  suggests  that  these  waves  originate  at  distances 
beyond  30  kilometers.  If  the  observed  turbulent  wake  results 
from  an  intrinsic  tendency  of  these  waves  to  shed  semiperiodic 
disturbances  behind  the  main  body,  or,  if  these  waves  are  breaking, 
then  they  are  rapidly  losing  energy  and  it  must  therefore  be 
anticipated  that  they  are  of  limited  range.  On  the  other  hand, 
if  the  observed  turbulence  results  from  the  disruption  of  a 
conditionally  unstable  atmosphere,  then  these  waves  could  have 
a range  of  up  to  about  600  kilometers,  the  limit  set  by  the 
speed  of  the  wave  and  the  duration  of  the  nocturnal  inversion. 

It  is  also  conceivable  that  these  waves  are  created  by  a different 
form  of  propagating  disturbance  which  interacts  with  the 
nocturnal  inversion. 
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Smart  (1966)  and  Jordan  (1972)  associated  their 
isur  l exponential  pr  ssure  pulses  wl  h thunderstorms.  The 
organized,  precipitation  initiated  and  maintained  downdraft  of 
cold  air  in  severe  mature  storms  is  well  documented  (Wallington, 
19ol ; Browning  and  Ludlam,  1962;  Spillane  and  McCarthy,  1969). 
Smart  developed  a model  to  describe  these  pressure  pulses  by 
considering  the  unstable  solution  of  Brunt's  equation  of  motion 
(Brunt,  1927)  for  the  vertical  displacement  (dh)  of  a parcel  of 


d2 (dh) 


|(g+  ^)dh  = 0 


(3.9) 


Brunt  noted  that  if  the  lapse  rate,  -dT/dh,  is  greater  than  the 
adiabatic  lapse  rate,  6 , no  restoring  force  exists,  and 
consequently  a downwards  displaced  air  parcel  will  continue 
to  accelerate  exponentially.  Smart  adapted  this  result  to 
explain  the  observed  form  of  the  pressure  pulses  by  noting  that 
the  downwards  descending  parcel  is  subject  to  a constant 
pressure  gradient  and  thus  the  density  of  the  parcel  increases 
exponentially  with  the  result  that  (assuming  the  approximate 
validity  of  the  hydrostatic  law)  the  ground  level  pressure 
contains  an  exponentially  increasing  component.  Smart  also 
hypothesized  that  the  downwards  motion  continues  to  accelerate 
in  smooth  laminar  flow  until  a critical  velocity  is  reached  at 
which  point  the  flow  becomes  turbulent  and  mixes  rapidly  with 
nearby  warmer  air  thus  causing  rapid  oscillations  behind  the 
pulse  and,  eventually,  a reduction  of  the  pressure  to  the  ambient 
level.  On  the  basis  of  this  model,  Smart  was  able  to  account 
for  4 out  of  10  observations  of  exponential  pressure  pulses. 
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March  are  due  to  thunderstorm  activity.  However,  only  32% 
of  all  events  that  have  been  observed  occurred  during  periods 
of  local  storm  activity.  A further  complication  arises  from 
the  observation  that  intense  local  storms  do  not  necessarily 
generate  solitary  atmospheric  waves.  For  example,  consider 
the  results  of  observations  made  during  the  month  of  January, 
a month  characterized  by  the  maximum  rainfall,  115-^  nun  over  a 
10  day  interval,  and  the  highest  number  of  counts  (1^90)  on  a 
lightning  flash  counter,  as  measured  at  the  Varramunga  Seismic 
Station.  During  this  period  only  one  very  weak  solitary  wave 
v/ as  recorded  and  this  particular  event  occurred  on  a day  when 
storm  activity  was  entirely  absent  throughout  the  Northern 
Territory. 

It  is  clearly  necessary  to  explore  other  possible  source 
mechanisms  for  the  solitary  atmospheric  w aves  at  VJarramunga. 

The  fact  that  these  waves  come  from  the  preferred  directions 
of  20°  and  1-0°  suggests  that  they  are  orographic  in  origin. 
However,  unless  one  is  willing  to  accept  the  hypothesis  that 
these  waves  originate  in  the  atmospheric  circulation  associated 
with  high  mountains,  specifically  the  high  mountains  of  New 
Guinea  and  the  Australian  Alps  (or  even  the  Southern  Alps  of 
New  Zealand)  - an  unlikely  hypothesis  for  which  there  is  no 
direct  evidence  - then,  the  fact  that  the  topography,  on  a more 
reasonable  scale  of  distance,  differs  considerably  in  these  two 
source  directions,  as  is  described  below,  indicates  that  two 
additional  mechanisms  may  be  involved  in  the  production  of 
solitary  atmospheric  waves  in  central  Australia. 

Consider  first  the  origin  of  isolated  waves  which  arrive 
from  an  aximuth  of  140°. 


These  waves  are  most  frequently 
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detected  near  midnight  but  they  may  also  occur  at  any  other 
time  during  the  night.  As  one  proceeds  in  the  direction  of 
140°  from  the  infrasonic  array  the  terrain  can  be  described  as 
a more  or  less  flat  featureless  plain  for  the  first  100 
kilometers.  The  Murchison  Range  begins  to  appear  5 kilometers 
to  the  south  of  this  path  at  about  10  kilometers  from  the 
array.  This  chain  of  low  hills,  which  runs  parallel  to  the 
140°  direction,  continues  to  rise  until  a maximum  elevation 
of  about  200  meters  over  the  surrounding  plain  is  attained  at 
a distance  of  about  80  kilometers.  Beyond  100  kilometers  the 
level  plain  slowly  gives  way  to  the  Davenport  Range  which 
intersects  the  Murchison  Range  at  the  same  altitude,  at  about 
140  kilometers  from  the  array.  Further  afield  in  this  direction 
the  terrain  descends  to  an  arid  plain  which  continues  uninter- 
rupted and  finally  merges  with  the  Simpson  Desert  at  a distance 
of  about  4-50  kilometers. 

■The  presence  of  these  nearby  ranges  suggests  that  the 
origin  of  solitary  waves  coming  from  this  direction  is  to  be 
found  in  the  interaction  with  an  existing  inversion  of  katabatic 
(down-slope)  nocturnal  flow  of  surface-cooled  air  from  these 
hills  on  to  the  surrounding  plain.  The  hypothesis  of  super- 
critical katabatic  drainage  which  leads  to  an  internal  hydraulic 
jump  has  been  employed  with  considerable  success  by  Ball  (1956) 
(see  also.  Ball,  1957;  Lied,  1964)  to  explain  the  sudden 
stationary  discontinuity  found  in  surface  winds  in  the  Antarctic. 
Similarly,  Clarke  (1972)  has  explained  a frequently  occurring, 
sudden,  near-dawn  squall  on  the  south  coast  of  the  Gulf  of 
Carpentaria  by  evoking  the  concept  of  a propagating  undular 
hydraulic  jump  created  in  the  process  of  nocturnal  surface 
flow  down  a 1:1000  slope  from  the  500  meter  highlands  to  the 


oast.  Clarke  carried  out  a number  of  numerical  experiments 
and  concluded  that  internal  atmospheric  hydraulic  jumps  of 
this  type  should  commonly  occur  at  low  latitudes  when  katabatic 
flow  discharges  on  to  a plain  under  conditions  of  a stable 
radiation  inversion. 

This  conclusion  is  reinforced  by  the  acoustic  soundings 
at  Julian  Creek  reported  by  Reynolds  and  Gething  (1970)  which 
revealed,  as  has  already  been  noted,  the  occurrence  of  a jump 
in  the  height  of  the  nocturnal  inversions.  These  authors  also 
hypothesised  the  existence  of  katabatic  flow  aided  in  part  by 
the  geotrophic  wind  and  orographic  convergence  and  concluded 
that  the  observed  jumps  in  the  inversion  level  could  be  attributed 
to  either  the  head  of  a gravity  current  propagating  as  an 
internal  bore  or  a travelling  internal  hydraulic  jump. 

It  is  worth  noting  at  this  point  that  none  of  the  data 
collected  so  far  from  the  infrasonic  array  indicates  the 
occurrence  of  either  a propagating  hydraulic  jump  or  the  frontal 
zone  of  an  advancing  atmospheric  density  current.  This  is  not 
an  unexpected  result  since  the  local  topography  indicates  that 
katabatic  currents  are  unlikely  to  flow  in  the  direction  of  the 
array.  It  is  therefore  proposed  that  the  observed  solitary 
v/aves  coming  from  an  azimuth  of  140°  are  generated  in  the 
interaction  of  katabatic  flow  with  the  noctural  inversion  and  that 
they  propagate  along  the  inversion  to  the  array.  It  is 
premature,  without  further  experimental  work,  to  examine  this 
mechanism  in  detail  but  it  probably  takes  one  of  the  following 
forms : 

(a)  the  direct  implusive  interaction  of  an  advancing 

gravity  current  impinging  on  a established  inversion; 


(b)  the  creation  of  a hydraulic  jump  at  the  point  where 
supercritical  down-slope  flow  encounters  the  horizontal 
plane  - this  hydraulic  jump  may  subsequently  propagate 
upstream  as  flow  conditions  vary; 

(c)  the  momentary  retardation  of  an  advancing  bore  upon 
encountering  a topographical  obstacle  such  as  a ridge 
or  valley,  and 

(d)  the  formation  of  a stationary  hydraulic  jump  in  conjunction 
with  a supercritical-subcritical  flow  transition  in  the 
neighbourhood  of  a topographical  barrier. 

The  interpretation  in  terms  of  katabatic  flow  at  an 
orographic  feature  accounts  for  both  the  occurrence  of  waves 
coming  from  the  direction  of  140°  and  for  the  scarcity  of  events 
originating  in  the  north  and  south-west  quadrants  where 
orographic  features  suitable  for  the  production  of  katabatic  flow 
are  absent.  It  is  possible  that  a few  events  originate  at 
long  range  in  the  Macdonnell  Ranges  to  the  south  and  south-west. 

An  attempt  to  apply  an  Interactive  katabatic  wind  inter- 
pretation to  account  for  the  source  of  the  large  number  of 
solitary  events  arriving  from  an  azimuth  of  20°  is  less 
successful.  In  this  direction,  the  elevation  decreases 
steadily  from  about  410  meters  MSL  to  200  meters  at  a distance 
of  about  150  kilometers  from  the  array.  From  this  point,  the 
land  rises  slowly  to  form  the  Barkly  Tablelands  at  an  altitude 
of  about  260  meters.  At  the  edge  of  the  Tablelands,  about 
340  kilometers  from  the  array,  the  land  again  rises  to  about 
310  meters  and  then  falls  continuously  to  the  hot  tropical  coast 
of  the  Gulf  of  Carpentaria,  500  kilometers  from  the  array. 
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It;  is  possible  that  solitary  waves  originating  in  this 
direction  are  produce!  In  the  interaction  of  katabatic  flows 
to  the  east  from  either  the  IloDougall  Ranges  on  the  outskirts 
of  Tennant  Creek  or  from  the  Short  and  Whittington  Ranges 
located  a further  50  kilometers  to  the  north.  This  inter- 
pretation is  unlikely,  however,  since  these  series  of  ridges 
rise  to  a maximum  of  only  50  meters  and  the  lateral  extent  of 
the  down-slope  run  to  the  surrounding  plain  would  seem  to  be 
too  limited  to  support  significant  katabatic  flow.  A different 
type  of  source  mechanism  is  therefore  required  to  explain  the 
solitary  waves  which  originate  in  this  direction. 

An  indication  of  the  range  of  isolated  waves  coning  from 
an  azimuth  of  20°  is  given  by  the  observation  that,  if  one 
ignores  all  storm-related  signals,  all  waves  of  this  type  arrive 
at  the  infrasonic  array  after  23.30  C.S.T.,  and  further,  that 
the  majority  of  these  events  are  observed  to  occur  around  03-00. 
This  suggests  that  the  source  of  these  waves  is  to  be  found  at 
a considerable  distance  from  the  array.  If  It  is  assumed  that 
these  waves  are  created  during  or  shortly  after  the  formation 
of  a stable  inversion  then  the  average  measured  speed  of 
10  meters/second  is  consistent  with  a source  in  the  neighbourhood 
of  the  coastal  edge  of  the  Barkly  Tablelands. 

It  is  therefore  tentatively  proposed  that  these  solitary 
waves  are  generated  through  the  interaction  of  a warm  upslope 
sea  breeze  circulation  from  the  Gulf  of  Carpentaria  with  a 

developing  cold  nocturnal  inversion  behind  the  rim  of  the 
Tablelands.  This  interaction,  if  not  direct,  might  take  the 
form  of  an  impulsive  addition  of  momentum  to  the  inversion  layer 
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during  the  collapse  of  stationary  orographic  lee  waves  or  during 
•he  dissipation  of  standing  eddies  following  the  late  afternoon 
transition  from  off-shore  flow  to  nocturnal  katabatic  drainage. 

This  concludes  the  description  and  interpretation  of 
atmospheric  internal  solitary  waves  which  are  observed  in  the 
arid  interior  of  Australia. 


» 


An  examination  of  the  properties  of  a total  of  85  waves  of  this 
type,  observed  over  a 10  month  period,  has  shown  that  these 
solitary  waves  propagate  along  the  nocturnal  inversion,  that 
they  arrive  at  the  infrasonic  array  from  two  different 
azimuths,  and  that  they  play  an  important  role  in  the  production 
of  atmospheric  turbulence.  It  is  concluded  that  the  solitary 
waves  observed  at  Warramunga  belong  to  a new  class  of  internal 
solitary  wave  discovered  by  Benjamin  and  by  Davis  and  Acrivos. 

A number  of  source  mechanism  which  may  lead  to  the 
production  of  solitary  waves  in  central  Australia  have  been 
proposed.  Nonetheless,  until  the  completion  of  further 
experimental  and  theoretical  work,  these  ideas  must  be  viewed 
as  conjecture.  It  is  hoped  that  an  expansion  of  the  Warramunga 
infrasonic  array  to  an  aperature  of  100  kilometers  will  provide 
the  necessary  data  required  to  determine  the  origin  of  these 
unusual  waves. 
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1.  Configuration  of  the  five-component  nicrobarograph 
army . 

2.  Measured  nicrobarograph  amplitude  response. 

3.  Internal  solitary  wave  surface  pressure  profiles 
corresponding  to  a 5°K  temperature  inversion  at  the 
950  mbar  level.  The  dimensionless  amplitude,  a , 
is  chosen  to  be  0.5- 

(a)  Classical  internal  solitary  wave  described  by  (3.5) 

(b)  Internal  solitary  wave  described  by  Benjamin's 
theory  (3-3)  far  a fluid  of  great  depth. 

i.  General  internal  solitary  wave  profile  shape  as  a function 

of  she  dimensionless  amplitude  a . Both  families  of  curves, 
are  plotted  in  the  dimensionless  quantities  X/H  and  Y/H 
where  H is  the  undisturbed  depth  of  the  lower  fluid  at 
infinity . 

(a)  Shallow  fluid  theory 

(b)  Deep  fluid  theory. 

5.  Examples  of  solitary  wave  spectra  recorded  at  the  Warramunga 
Seismic  Station.  These  well-defined  waves  of  elevation 
are  typical  of  most  of  the  large  amplitude  events  which  have 
been  observed. 

6.  Further  examples  of  solitary  waves  of  elevation.  Note  the 
sudden  change  in  the  character  of  the  micropressure  spectrum 


following  the  passage  of  the  solitary  wave  illustrated  In  (a). 


7-  Various  forms  of  solitary  wave  spectra.  Figure  (a)  in  this 
group  is  chosen  to  illustrate-  an  example  of  a solitary  wave 
of  elevation  followed  by  a solitary  wave  of  depression. 

8.  Examples  of  solitary  waves  of  elevation  which  preceed  a 

decaying  train  of  coherent  long  period  spectral  components. 

9*  Examples  of  solitary  waves  of  elevation  which  mark  the  onset 
of  long  period  incoherent  components  in  the  atmospheric 
spectrum. 

10.  Examples  of  microfc arograph  records  which  illustrate  a 
variety  of  unusual  solitary  wave  profiles. 

11.  Frequency  of  observation  of  solitary  wave  signals  as  a 
function  of  calendar  day. 

12.  Frequency  of  observation  of  solitary  -waves  as  a function 
of  time  of  day. 

13.  Occurrence  of  solitary  waves  as  a function  of  source 
azimuth. 

It.  Diurnal  variation  of  the  micropressure  spectrum  illustrating 
the  growth  and  dissipation  of  the  nocturnal  inversion  and 
the  influence  of  solitary  wave  activity. 

15.  Some  examples  of  the  influence  of  solitary  waves  on  the 
atmospheric  spectrum. 

16.  A comparison  of  solitary  wave  pressure  profiles  for  the  case 

of  (a)  a wave  with  a relatively  quiet  wake,  and  (b)  a wave  which 
preceeds  a period  of  significant  atmospheric  turbulence. 
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